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ABSTRACT 

Cosmic microwave background (CMB) radiation data obtained by different experi- 
ments contain, besides the desired signal, a superposition of microwave sky contribu- 
tions. We present a fast and robust method, using a wavelet decomposition on the 
sphere, to recover the CMB signal from microwave maps. An application to WMAP 
polarization data is presented, showing its good performance particularly in very pol- 
luted regions of the sky. The applied wavelet has the advantages of requiring little 
computational time in its calculations, being adapted to the HEALPix pixelization 
scheme, and offering the possibility of multi-resolution analysis. The decomposition 
is implemented as part of a fully internal template fitting method, minimizing the 
variance of the resulting map at each scale. Using a characterization of the noise, 
we find that the residuals of the cleaned maps are compatible with those expected 
from the instrumental noise. The maps are also comparable to those obtained from 
the WMAP team, but in our case we do not make use of external data sets. In addi- 
tion, at low resolution, our cleaned maps present a lower level of noise. The E-mode 
power spectrum Cf^ is computed at high and low resolution; and a cross power spec- 
trum Cj^ is also calculated from the foreground reduced maps of temperature given 
by WMAP and our cleaned maps of polarization at high resolution. These spectra 
are consistent with the power spectra supplied by the WMAP team. We detect the 
E-modc acoustic peak at ^ 400, as predicted by the standard AC DAI model. The 
B-mode power spectrum Cf^ is compatible with zero. 

Key words: methods: data analysis - cosmic microwave background 



1 INTRODUCTION 

Component separation is a critical aspect in the analy- 
sis of cosmic microwave background (CMB) data. A good 
characterization of the data is a prerequisite to the ade- 
quate estimation of cosmological parameters. This need be- 
comes crucial when, as happens in B-mode detection ex- 
periment s, foreground a mplitudes are well above the sig- 
nal (e.g.. iTucci M. et al.| [2005V Two physical galactic pro- 
cesses are the major contaminants to CMB polarized sig- 
nal: synchrotron radiation and thermal dust. Both appear 
at large scales, are highly anisotropic and the spatial vari- 
ation of their emissivity is smooth. Besides, extragalactic 
emission also contaminates this cosmological signal: point 
sources and clusters are compact objects, roughly isotropi- 
cally distributed in the sky and every single object has a par- 



e-mail: cobos@ifca.unican.es 



ticular frequency dependence. Most of the component sepa- 
ration methods take into account only diffuse components, 
assuming that we are previously masking the brightest point 
sour ces or subtracting them by, typically, fitting approaches 
(seell erranz fc Vielvallioiol . for a recent review). 



C urrent and future e xperiments llRubifio-Martm et al 



20081 : iBrown et al] l2009l : ISievers et aLi[2009l: [Arnold et al 



l2006l : ICrainger etsl\ |2008| : 



20 la iKogu t et al.' 
Charlassier et al. 2008|) are able to measure CMB po- 
larization anisotropics with such precision that foreground 
contamination have become the major limitation when we 
try to analyze the data. This is the principal reason to invest 
effort and time in developing new techniques for separating 
components. The goal of all the proposed methods is to 
separate or, at least, to identify CMB anisotropics from the 
other components. The range of proposals includes internal 
linear combinations (ILC), Bayesian methods and indepen- 
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dent component analysis (see lDelabrouille fc Cardosdl200'7l , 
for a recent review). 

There is abundant literature that includes ap- 
plications of the various methods related to some 
polarization experiment s in vogue as, for instance, 



PLANCK (Leach etal 



iBetoule et al 



llGold et al.l 



vogu 

■ .20081: lEfstathiou et ahl [2009; 
20091 : iBaccigal upi et all 120041') a r id WMAP_ 
201ll : [pclabrou ille J. et al.1 |2009|: iKim et al ' 

Hi 



2007: Eriksen et al 



2009; Bona ldi et all 120071 : IMaino et al, 
200& .,2008). 

The method that we present in this paper is situated 
in the context of the internal linear combinations and it is 
a bet for a template cleaning in which coefficients are fitted 
in the space of a particular wavelet that enables a multi- 
resolution analysis. A fitting by scales allows, in practice, 
some efi'ective variation of the coefficients in the sky, which is 
an advantage over the template cleaning in real space. In this 
sense, our approach based on wavelet space effectively lies in 
between standard linear combination techniques applied in 
the r eal space and more sophi s ticated parametric methods 
fe.g.. lEriksen et al.ll2006l . l2008l : IStivofi et al.l[201ol l. 

Our approach is a fast procedure that especially shows 
its effectiveness in polluted regions, such as those that ap- 
pear in polarization experiments. 

This paper is structured as follows. The methodology 
is described in detail in Section [2l We set out an analysis of 
the low-resolution polarization WMAP data in Section^ In 
section[4l we show the treatment for high-resolution WMAP 
data in order to obtain the Cf^ and Cj^ spectra. Finally, 
we present the conclusions and discussion in section [S] 



2 METHODOLOGY 

In this work, we present a multi-resolution internal template 
cleaning (MITC) method for foreground removal. This is 
the initi al step of the map cleaning pro c ess in the SEVEM 
method (|Martmez-Gonzalez et~al]|2003l : iLeach et al.1 120081 ') 
to the case of polarization. 

For many purposes, it is a key point to have CMB maps 
at several frequencies instead of a single map. For instance, 
it would serve as a consistency check to verify whether any 
detected feature of the data is actually monochromatic or 
not (as, for instance, the case for non-gaussianity analysis). 

Another advantage of the method is that we do not 
need a thorough knowledge of foregrounds, because we ob- 
tain all the information to construct different templates 
from the data. Furthermore, this procedure preserves the 
original resolution of the CMB component. But the down- 
side is that the internal templates are noisy, so we in- 
crease the total noise level when we remove them from 
the data. This circumstance results, for instance, in an in- 
crease in the error bars of the power spectrum at high mul- 
tipoles. An alternative would be to incorporate external 
templates, created from data from other independent ob- 
servations or based on theoretical arguments. However, the 
current knowledge of foreground emissions, in polarization, 
is not substantiated with suitable ancillary data set, and 
for that reason, this option is not considered in this case. 
This situation may change in the future with the informa- 
tion expected to be provided by experiments like PLANCK 
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Figure 1. Outline of construction of the detail coefficients at 
resolution j (di) as the substraction of the approximation coeffi- 
cients at resolution j-1 (yi), from the approximation coefficients 
at resolution j (xi). 



llTauber et al."20 10l). C-BASS (iKing et al.l l201ll) or QUI- 
JOTE l(Rubifio-Martm et al.ll2008t ). 



2.1 The HEALPix wavelet 

Wavelets are a powerful tool in signal analysis and 
are extensively used in many astrophysics applications. 
Several examples of implementation of component sep- 
aration methods which employ ver y diverse wavelets 
can be found in th e liter at ure (e.g.. [ Ghosh et al] 
Delabrouille J. et al.l l2009l: iGonzalez-Nuevo et al 



2011 



2006 



Vielvaet al.ll2003l : [Hansen et al.l 120061 ). Thev are localized 
wave functions, that allow for a multi-resolution treatment 
of the data. This fact represents an advantage over other 
component separation methods because it allows us to vary 
the effective emissivity of foregrounds. 

We use the s o-call ed HEALPix wavelet, HW, 
ijCasaponsa B. et al.l l201ll ). a discrete and orthogonal 
wavelet that provides a multi-scale decompositio n on the 
spher e adapted to the HEALPix pixelization ( | G6rski et al.l 
l2005h . The resolution of a map in the HEALPix tessellation 
is given in terms of the Ngide parameter, defined so that the 
number of pixels needed to cover the sphere is A*' = 12A^j^g . 
The resolution j of a map is a number such that 2^ = Nstde- 
A CMB map is decomposed in the wavelet coefBcient space 
in a series of maps from the resolution of the original map 
to the lowest resolution considered. All of these maps, ex- 
cept the lowest resolution one, are called details. The last 
one is called the approximation, and is constructed by de- 
grading the original map to the appropiate resolution, i.e., 
to calculate the approximation coefficient at resolution j-1 
at a given position i we take the average of the four daugh- 
ter pixels at resolution j. The way that different detail maps 
are built is illustrated in figure [1] At each resolution j, detail 
coefficients are calculated as the substraction of the approxi- 
mation coefficients at resolution j-1 from the approximation 
coefficients at resolution j. Both this process and the math- 
ematical formalism o f this wavelet is carefully explained in 
ICasaponsa B. et al.l (|201ll ). where the HW is used to put 
constraints on the /nl parameter from WMAP data. 

In that paper, it is said that the reconstruction of a map 
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M{xi) can be written as 

p=0 

4 J-1 
m=l i=jo P=0 



(1) 



where Ajp 



and 'ym,j,e are the approximation and detail co- 
efficients respectively, 4ij,p{xi) is the scaling function and 
ipm,j,p{xi) refers to the wavelet functions. The j index takes 
values from the highest resolution J to the approximation 
resolution jo. 

The advantage of this wavelet with respect to others, in 
addition to its straightforward implementation, lies in the 
speed of the involved operations. The computational time 
for the wavelet decomposition is of the order of the num- 
ber of pixels (~ Npix) whereas, for example, for the conti- 
nous wavelet transform of the sp herical Mexican hat wavelet 
iMartmez-Gonzalez et alll2002l ) or the needlets (|Baldi et al.l 
I2OO9I ') this time is of the order of - N. 



3/2 



2.2 Template fitting 

The signal Tj{p) at resolution j is constructed by subtract- 
ing a linear combination of different templates tij from the 
original signal, Tj, as follows 



Nt 



dp), 



(2) 



where Nt is the total number of templates and p is a pixel 
index. 

An internal template is formed as the difference of two 
maps of the same resolution, corresponding to different fre- 
quencies, in units of thermodynamic temperature. 

The variance of the cleaned map is optimally minimized 
at each scale to obtain the coefficients Pij or, equivalently, 
the quadratic quantity 



Xj = 



^{f,(p)C-^[f,(p)]*} 



(3) 



where is the inverse of the covariance matrix calculated 
as the sum of contributions of the CMB and instrumental 
noise (both, from the map to be cleaned and the templates). 

From the previous discussion, it is obvious that the ap- 
proach to produce an optimal recovery of the CMB would 
require a certain knowledge of this signal, via its correla- 
tions. However, a more robust estimator, without a priori 
knowledge of the signal to be estimated, may be built by 
considering only the instrumental noise correlations. 

We have checked, however, that the gain in the CMB 
recovery, by including the information related to the instru- 
mental characteristics is, in practice, very little. Even more, 
in some situations (as it is the case of the WMAP full resolu- 
tion data, see section [4| the instrumental noise information 
is limited to the autocorrelation. Therefore, in this work we 
have decided to perform the internal template fitting with 
uniform weights for all the pixels at each scale, which implies 
to minimize the following quantity: 



E 



Nt 



(4) 



Finally, we recover a single map performing the wavelet 
synthesis. It can be written as 



f{x) = T(£) -EE 7..(^)i..(^), 
i=i j=i 



(5) 



where Nres denotes the number of involved resolutions and 
7ij some new coefficients given as linear combinations of Pij 
coefficients which are the result of the synthesis process. 



3 ANALYSIS OF LOW RESOLUTION WMAP 
DATA 

The instrument al noise in WMAP polarization is known 
to be correlated (|jarosik et al.ll201ll '). Ahhough the WMAP 
data are typically given at a HEALPix resolution of N^ide = 
512, a more accurate version of the pixel-to-pixel correla- 
tion is only available at low resolution, namely, A^side = 16. 
Taking into account this difference, we have performed the 
cleaning of the WMAP data in two cases: for low and high 
resolution maps. In this section, we analyse the maps at 
Nside = 16. 

The WMAP data are composed by, at least, a superpo- 
sition of CMB, synchrotron and thermal dust emissions. The 
WMAP team proposed a template fitting in the pixel space 
to clean the foreground emission in the Ka, Q, V and W 
maps, using as templates the K ba nd (for the synch r otron ) 
and a low resolution version of the iFinkbeiner et all (|l999l ) 
model for the thermal dust, with p olarization direct ion de- 
rived from starlight measurements (|Gold et al.ll20lj ). 

In our approach, we use only a synchrotron template, 
constructed as K-Ka. The reason for neglecting the thermal 
dust template is because a previous analysis in real space 
shows that its coefficients are much smaller than the cor- 
responding ones for the synchrotron template. We clean Q 
and U polarization components independently minimizing 
the variance of the cleaned maps of the Ql, Q2, VI and V2 
differencing assemblies (DAs). The wavelet decomposition 
is carried out down to resolution j = 3 for the data map, 
thus, in addition to the approximation, we have a single de- 
tail map at j = 4. Best fitting coefficients for the considered 
DAs are given in table[T] We apply the WMAP polarization 
analysis mask that excludes a 26% of the sky. 



3.1 Cleaned maps 

Since the CMB polarization signal is clearly subdominant 
in the WMAP low resolution data, it is hard to establish a 
criterion to evaluate the goodness of the cleaning process, 
and to perform comparisons with different solutions. 

We have decided to evaluate this goodness by compar- 
ing the cleaned map with the expected signal for a noisy sky 
following the WMAP instrumental noise characteristics. In 
this sense, a good compatibility with the noise properties 
would indicate that foregrounds have been satisfactorily re- 
duced. 

We generate a set of 10* simulations of the noise maps 
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Frequency band 




Ql 


Q2 


VI 


V2 


Detail (j = 4) 


Q Stokes 
U Stokes 


0.092 
0.074 


0.103 
0.117 


0.036 
0.020 


0.023 
0.048 


Approximation (j = 3) 


Q Stokes 
U Stokes 


0.244 
0.241 


0.259 
0.236 


0.081 
0.085 


0.125 
0.112 



Table 1. Template cleaning coefficients for Q and U Stokes parameters and DAs for the low resolution case. 




Figure 2. Upper panels show the distributions of our cleaned maps and the bottom ones present the same results for the WMAP 
procedure. The solid line (red) corresponds to the theoretical curve of a with as many degrees of freedom as pixels outside the mask 
(i.e., the effective number of pixels in Q and U maps: 4518). The dashed line is the distribution calculated from simulations of our cleaned 
maps. The vertical line shows the value of the data maps in each case. The columns correspond to different frequency bands, from 
left to right: Q, V and W bands. 



resulting from our MITC method at Q, V and W frequency 
bands, Mr{p), with r £ {1, 10''}, in order to construct a 
distribution, calculating each value as 

xl^Y. Mr{p)NoUp, q)MUq), (6) 

where Nobs is the noise correlation matrix. A number of 
simulations of the order of a million is required to estimate 
this matrix so that the distribution converges to the the- 
orical curve of a distribution with as many degrees of 
freedom as pixels outside the mask in Q and U maps (in 
this case, we have 4518 degrees of freedom). This distribu- 
tion characterizes the expected noise level at each frequency 
map. We can associate the value of the data with relative 
levels of signal. We can say that the cleaned maps contain 
more than just noise (typically foreground residuals, since 
the CMB is subdominant compared to the noise at these 
scales) if the data value is much higher than typical values 
of the distribution. Conversely, we can ensure that our maps 
are compatible with the expected noise and that residuals 
are small if the data value falls within the distribution. The 



Frequency band 


Q 


V 


W 


^ Forered 


4566 
4709 


4453 
4586 


4762 
4787 



Table 2. Different values of computed with our seven- 
year cleaned maps per frequency band iXcieaned^ ^'^'^ with 
WMAP seven-year foreground reduced maps per frequency band 

orered) ■ 



DA 


Ql 


Q2 


VI 


V2 




4489 


4546 


4405 


4486 



Table 3. Different values of x^ computed with our seven-year 
cleaned maps per DA. 



values for each band are listed in table [2] and for each 
DA in table 131 

Our test is based on the assumption that the CMB con- 
tribution is negligible. We have tested that the CMB pro- 
vides a very small contribution (a shift of ~ 10 units of x^) 



(xi 
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-3 

X 10 

W2-W4 




2 

X 

Figure 3. The combinations of the raw W-band maps. Theo- 
retical curve of the distribution is represented by a solid red 
line and the values are shown by successive vertical lines, from 
left to right: W2-W4, W1-W2, (W1+W3)-{W2+W4), W2-W3, 
W1-W4, W3-W4 and W1-W3. 

by generating 10'* simulations with CMB and instrumen- 
tal noise of the cleaned maps. These simulations have been 
used to compute another distribution with the matrix 
that we have already calculated with only the noise compo- 
nent. When distributions are compared with each other we 
observe this typical deviation. Thus, the CMB contribution 
to the value of the of the data is negligible and, there- 
fore, any significant deviation from the mean value has to 
be assigned to foreground residuals. 

An indirect comparison can be made between the 
WMAP procedure and our MITC method through the rela- 
tive possitions of the value of the data with respect to the 
distribution. As seen in figure [H we obtain that the value 
of the cleaned maps is fully compatible with instrumental 
noise at Q and V frequency bands. At W band the x^ value 
of the cleaned data is in the tail of the distribution probably 
due to the presence of foreground residuals. The deviation 
is even larger when the WMAP procedure is used. A sig- 
nificant improvement is also found at Q band since the x^ 
value is shifted from 2a to 0.5a when our MITC method is 
used. 

In addition, although we use a template that is noisier 
than the ones used by the WMAP team, the noise levels of 
our cleaned maps are lower. We have measured a difference 
of about a 10% in terms of the standard deviation of the 
data maps (this difference is confirmed by instrumental noise 
simulations) . 

In order to check further the apparent excess of signal 
at W band obtained by the two approaches, we computed 
analytically the noise covariance matrix of different com- 
binations of the raw W-band DAs maps which contain, in 
principle, only a combination of instrumental noise. With 
these covariance matrices, based on the full-sky covariance 
matrices of each DA, a x^ value of the data maps is obtained. 
It is shown in figure [3] that these maps are still compatible 
with the expected noise. However, it is significant that all 
values are to the left of the distribution and that the most 
deviated ones involve W2, followed by W4. We have also 



analysed the distribution of the x^ values of the single-year 
foreground reduced maps supplied by the WMAP team for 
each DA at the W band, obtaining values more deviated to- 
wards the tails for the W2 and W4 DAs. This may suggest a 
not good enough characterization of the instrumental noise 
for these DAs. 

3.2 Polarization power spectra, Cf^ and Cf^ 

We carry out an estimation of the polarization spectrum us- 
ing our cleaned maps of the Ql, Q2, VI and V2 DAs. A 
pseudo cross-power spectrum Df^ between any two differ- 
encing assemblies A and B can be calculated as 

Df^ = Y: Mff \p, \'B^Bf, {Cr) + {Nn, (7) 
t' 

where A,B & {Ql, Q2, VI, V2\ A^ B}; and, in the case of 
an EE power spectra, 

riAB _ ^ \ ^ A B* fQ\ 

where eirn are the spherical harmonic coefficients of the E- 
mode. Assuming a circular beam response, we denote the 
beam of the A map as Bf and the window function of the 
HEALPix pixel by pe; (iV/^> is the noise cross-power spec- 
trum. The bias introduced by this term comes from the inter- 
nal template fitting procedure. It is small and controled by 
simulatio ns. Finally, the cou pling kernel matrix M^^i is de- 
scribed in lHivon et al.l (|2002D and, for the case of th e polar - 
ization components, in Appendix A of iKogut et al.l (|2003h . 
This procedure is usually referred to as MASTER estima- 
tion. An estimator, Ce, can be computed as a linear combi- 
nation of the six different spectra weighted by the inverse of 
their variances in the following way: 

where i = AB and af — aAas- These variances are given 
by the WMAP team in the LAMBDA web siteQ- 

The resulting power spectra are shown in figure]?] From 
the Cf^ spectrum we can say that most of the values are 
compatible with zero, so there is almost no signal except 
perhaps for low multipoles £ < 6. As expected, the B-mode 
spectrum Cf^ signal is compatible with zero. Both spectra 
are compatible with those that the WMAP team supplies. 
Our error bars are larger than those obtained by the WMAP 
team, because of the use of an estimator that is not optimal, 
a pseudo-spectrum, whereas the WMAP team uses a pixel- 
base likelihood. 



4 ANALYSIS OF HIGH RESOLUTION WMAP 
DATA 

In this section, we analyze WMAP data maps at Nside = 
512. This approach allows us to study the cleaning at smaller 
scales where, a priori, the correlation of the noise is less im- 
portant. So then we only take into account the noise covari- 
ance matrix of each pixel. In this case, the cleaning method 

^ http://lambda.gsfc.nasa.gov/ 
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Figure 4. Polarization power spectrum EE (upper pannel) and 
BB (bottom pannel) for low resolution analysis. Circles are the 
spectrum supplied by WMAP team and asterisks represent our 
estimation. The fiducial model is plotted by the solid line. 



based on the wavelet space is applied using two different in- 
ternal templates. The first one is constructed as K-Ka and 
accounts for the synchrotron radiation. The second one is 
built as VI- W3 and attempts to characterize the thermal 
dust. The V and W DAs to be cleaned have been selected 
by having a lower noise. As in the low resolution case, the 
wavelet decomposition is carried out until resolution j = 3 
for the approximation map, hence we have additionally 6 
different detail maps in this high resolution case. Again, the 
WMAP polarization mask is used. Best fitting coefficients 
are shown in table ID 

From the cleaned Q and U maps we study the power 
spectra, Cf^ and Cf^ , as in the previous section. The spec- 
trum error bars are estimated from 10^ noise simulations. In 
general, the CMB contribution is neglegible compared to the 
noise one. Bins are taken as a weighted average of the mul- 
tipoles involved. These weights are calculated as the inverse 
of the variance of each Ce. 

The spectra are compatible with those that the WMAP 
team has obtained. Our error bars are larger at high multi- 
poles, since, at this scale range, the number of effective cross- 
spectra is much smaller than the one used by the WMAP 
team (where all the W-band DAs are available). Neverthe- 
less, as seen in figure [51 our cleaned maps seem to present 
a lower level of contamination than those supplied by the 
WMAP team. We expected a better foreground removal 
since the templates used are closer to the foreground signal 
distribution accross the sky in our case. However, whether 
this may have an impact on the determination of the cos- 
mological parameters is not clear and would require an ex- 
haustive analysis which is outside the scope of this paper. 

Finally, the two points accounting for the largest scales 
of the spectra are taken from figure U since the correlation 
of the noise at these scales is important, and it has been bet- 
ter modeled in the previous section, where a more accurate 
version of this information was available. 




30 - 
20 - 
10 - 
i 



-10 - 
-20 - 
-30 - 
-40 - 
-50 - 



1 

D 
-1 

-2 



Figure 5. Polarization power spectrum EE (upper pannel) and 
BB (bottom pannel) for high resolution analysis. The circles are 
the spectrum supplied by the WMAP team and the asterisks 
represent our estimation. The solid line represents the fiducial 
model for CF^ and the zero value for C?^ . 



The resulting power spectra are presented in figure [S] 
Our outcome is compatible with the WMAP team analysis, 
where it is possible to distinguish the acoustic peak around 
£ ~ 400 in the E-mode spectrum. As expected, the B-mode 
power spectrum is compatible with zero. 

Our independent approach can be seen as a confirma- 
tion of the previously detection reported by the WMAP 
team. 

A similar procedure is applied to determine the corre- 
lation between temperature and E-mode polarization data. 
We carry out the analysis using our cleaned maps of the 
Q and U Stokes parameters and the foreground reduced 
temperature maps that the WMAP team supplies. In this 
case, the combination of two maps of the same DA is al- 
lowed and the equation [7] has the same form as long as 
we add that, in Cf^, A refers to the DA of temperature 
maps and B to the DA of polarization maps, with A,B £ 
{(51,(32,^2,^1,14/2,14/4}. For the temperature maps, we 
use the temperature analysis mask that the WMAP team 
supplies a nd the MASTER es timation is computed as is de- 
scribed in lKogut et"aLl (120031 ). 
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DA 




Template 


Stokes 


Ql 


Q2 


V2 


Wl 


W2 


W4 


Detail {j 


= 9) 




Q Stokes 
U Stokes 


-0.0175 
0.0525 


-0.0067 
0.0173 


0.0024 
0.0085 


0.0754 
-0.0869 


0.0072 
0.1171 


-0.0222 
-0.0562 


(V1-W3) 


Q Stokes 
U Stokes 


-0.0008 
0.0007 


0.0023 
0.0006 


-0.0015 
0.0022 


0.0013 
0.0005 


-0.0015 
-0.0226 


0.0054 
-0.0106 


Detail (j 


= 8) 


(K-Ka) 


Q Stokes 
U Stokes 


-0.0238 
-0.0080 


0.0185 
-0.0176 


0.0135 
-0.0059 


-0.0046 
-0.0513 


0.0076 
0.0138 


0.0311 
-0.0196 


(V1-W3) 


Q Stokes 
U Stokes 


-0.0044 
0.0012 


-0.0003 
0.0028 


-0.0008 
0.0042 


-0.0048 
0.0013 


0.0020 
0.0030 


0.0071 
0.0012 



Detail (j = 7) 



(K-Ka) 
(V1-W3) 



Q Stokes 
U Stokes 
Q Stokes 
U Stokes 



0.0006 
-0.0208 
0.0054 
-0.0017 



0.0045 
0.0196 
0.0004 
0.0006 



-0.0091 
0.0038 
-0.0002 
0.0032 



-0.0227 
-0.0263 
-0.0053 
0.0019 



0.0288 
0.0047 
-0.0002 
-0.0036 



-0.0224 
-0.0142 
0.0069 
0.0009 



Detail (j = 6) 



(K-Ka) 
(V1-W3) 



Q Stokes 
U Stokes 
Q Stokes 
U Stokes 



0.0224 
0.0278 
0.0079 
0.0066 



-0.0112 
0.0172 
0.0014 
-0.0051 



-0.0073 
0.0104 
-0.0039 
0.0016 



0.0302 
0.0144 
-0.0070 
-0.0012 



-0.0057 
0.01797 
-0.0095 
-0.0096 



0.0407 
0.0076 
-0.0038 
0.0179 



Detail (j = 5) 



(K-Ka) 
(V1-W3) 



Q Stokes 
U Stokes 
Q Stokes 
U Stokes 



0.0351 
0.0567 
-0.0131 
-0.0080 



0.0729 
0.0523 
0.0012 
0.0009 



0.0702 
0.0670 
0.0207 
-0.0071 



-0.0325 
0.1166 
-0.0070 
-0.0362 



0.0253 
0.0255 
-0.0578 
-0.0006 



-0.1171 
-0.0225 
0.0508 
0.0363 



Detail (j = 4) 



(K-Ka) 
(V1-W3) 



Q Stokes 
U Stokes 
Q Stokes 
U Stokes 



0.2292 
0.0741 
-0.0708 
-0.0440 



0.0571 
0.1172 
0.0774 
0.0359 



0.0548 
-0.0844 
-0.1389 
-0.0397 



-0.0835 
0.1445 
0.0386 
0.1985 



0.3029 
-0.2241 
0.1316 
-0.2028 



-0.2239 
-0.2723 
-0.0135 
0.2320 



Approximation {j = 3) 



(K-Ka) 
(V1-W3) 



Q Stokes 
U Stokes 
Q Stokes 
U Stokes 



0.1732 
0.1397 
-0.2852 
0.2087 



0.2976 
0.2315 
0.2030 
0.1022 



0.1988 
0.1674 
0.6050 
-0.0048 



0.1413 
0.1607 
0.5513 
0.5052 



0.1234 
0.4529 
0.2601 
0.0160 



0.2612 
0.5386 
-1.1963 
-0.3867 



Table 4. Template cleaning coefficients for Q and U Stokes parameters and different frecuency bands for the high resolution case. 




1 2 

10 10 
1 



Figure 6. Polarization power spectrum TE. Circles are the spec- 
trum supplied by WMAP team and asterisks represent our esti- 
mation. The fiducial model is plotted by the solid line. 



Since the CMB cosmic variance in temperature con- 
tributes significantly, we cannot ignore it this time in the 
calculation of the error bars. So we generate 10^ simula- 
tions of CMB plus noise, which undergo the same process of 
cleaning and combination to get the cross power spectrum, 
with which we estimate the error bar. These errors are rela- 
tively larger with respect to the WMAP results, due to the 
less number of effective cross spectra. In addition, we re- 
mark that, while we use a pseudo-spectrum, low resolution 
analysis is performed through a pixel-base likelihood by the 
WMAP team. 



The resulting cross power spectrum is presented in fig- 
ure [S] and it is compatible with the power spectrum that 
the WMAP team supplies. 



5 CONCLUSIONS 

We introduce an internal template cleaning method that 
uses a wavelet decomposition on the sphere. Among its ad- 
vantages, it is included the possibility of multi-resolution 
analysis, allowing an effective variation of the spectral index 
in the sky. Much lower computational time is needed than 
with other widely-used continous wavelets. In addition, a 
good treatment of incomplete sky coverage is given because 
of the compact support of the HEALPix wavelet. 

The MITC method result is a set of some cleaned maps 
at several frequencies that can be used, for instance, to ver- 
ify whether any detected feature of the data is actually 
monochromatic or not. The exclusive use of internal tem- 
plates allows us to analyze the maps without making any 
prior assumptions about the foregrounds in polarization. 
However, although the implementation that is shown in this 
work make use of only internal templates, it can be trivially 
extended to deal with external templates as well. 

We perform an analysis of 7-year WMAP data obtain- 
ing outcomes that are compatible with the WMAP team 
results. Furthermore, we have hints of better cleaning of the 
Q-band map at large scales and we have obtained cleaned 
maps, at least, as good as those that the WMAP team sup- 
plies for V and W bands. Let us remark that our approach 
does not make use of any additional template: everything is 
obtained from WMAP data. We have checked that, although 
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we use noisier templates, the instrumental noise levels of the 
final cleaned maps are lower than those of the maps provided 
by the WMAP team. 

High resolution maps are also analysed. In agreement 
with the WMAP team, we find an E-mode detection at ^ ~ 
400, as predicted by the standard ACDM model. We also 
obtain that the B-mode level is compatible with zero. These 
independent findings are a confirmation of the result already 
presented by the WMAP team. 

The clean maps produced at this work, both at low 
and high resolution, are available at the following website: 
|http://maj{:.ifca.unican.es/cobos/WMAP 7yrPOL^ 
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